Infinite invariant densities for anomalous diffusion in optical lattices and other 

logarithmic potentials 
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We solve the Fokker-Planck equation for Brownian motion in a logarithmic potential. When 
the diffusion constant is below a critical value the solution approaches a non-normalizable scaling 
state, reminiscent of an infinite invariant density. With this non-normalizable density we obtain the 
phase diagram of anomalous diffusion for this important process. We briefly discuss the consequence 
for a range of physical systems including atoms in optical lattices and charges in vicinity of long 
polyelectrolytes. Our work explains in what sense the infinite invariant density and not Boltzmann's 
equilibnum describes the long time limit of these systems. 
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Brownian particles in contact with a thermal heat bath 
and in the presence of a binding potential field V{x) at- 
tain a steady state which is the Boltzniann equilibrium 
distribution Wsteady{x) = N exp[—V{x)/kbT]. An inter- 
esting case is the logarithmic potential: V{x) oc ep ln(a;) 
for a; 3> 1. Inserting V{x) into the Boltzmann distribu- 
tion, one finds that the steady state solution is described 
by an asymptotic power law Wsteady{x) ~ Nx~'^°^'''''^ . 
eo/kbT must be larger than one for the normalization N 
to exist. Brownian motion in a logarithmic potential has 
attracted much attention since it describes many physi- 
cal systems, ranging from diffusion of momentum of two 
level atoms in optical lattices [l|-|3|, singleparticle mod- 
els for long ranged interacting systems [J, |5| , the famous 
problem of Manning condensation describing a charged 
particle in the vicinity of a long and uniformly charged 
wire (i.e. an ion in vicinity of a long charged polymer) 
[6| , and very recently the motion of nano-particles in an 
appropriately constructed force field [7|. 



In this Letter we provide the long sought after |1|- 
long time solution of the Fokker-Planck equation describ- 
ing the dynamics of Brownian particles in a logarithmic 
potential. Naively one would expect that in the long 
time limit the equilibrium distribution, i.e. the Boltz- 
mann distribution, is reached; however the logarithmic 
potential turns out to be much more interesting. To start 
with, we point out that the second moment in the steady 
state {x'^)steady may diverge, namely {x'^)steady = oo if 
1 < eo/khT < 3. This behavior is unphysical. In par- 
ticular in the context of optical lattices, it implies that 
the averaged kinetic energy of the atoms is infinite which 
is of course impossible (see details below). If we view 
the problem of Brownian diffusion in a logarithmic po- 
tential dynamically, we immediately realize that the pro- 
cess cannot be faster than diffusion; namely (cc^) < 2Dt 
where D is the diffusion constant so (a;^) = oo is wrong 
(we will soon derive this bound from the Fokker-Planck 
equation). In this sense the steady state solution, e.g. 
the Boltzmann distribution, for a particle in a logarith- 
mic potential, does not describe the statistical properties 



of the problem, even in the limit of long times. 

Here we show that Brownian particles in a logarith- 
mic potential are characterized by an infinite invariant 
density. This density is not normalizable (hence the term 
infinite); however as we show, it does describe the anoma- 
lous behavior of the system. For example it can be used 
to obtain correctly the moments of the process, while the 
normalizable Boltzmann distribution completely fails to 
do so. We examine these issues first in the context of 
diffusion of momenta of atoms in an optical lattice, since 
this system is an excellent candidate to experimentally 
test our predictions in the lab. The reader should note 
that our results with some small notational changes de- 
scribe a wide class of Brownian trajectories in the pres- 
ence of a logarithmic potential (see discussion below). 

Fokker-Planck equation. The equation for the proba- 
bility density function (PDF) W{p, t) of the momentum 
p of an atom in an optical trap is modeled within the 
semi-classical approximation according to [l|-|3| 
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restores the momentum to its minimum while D describes 
stochastic momentum fluctuations which lead to heat- 
ing. From the Sisyphus effect, interaction of atoms with 
the counter propagating laser beams means that D is 
determined by the depth of the optical potential [l|43|, 
which in turn leads to experimental control of the un- 
usual statistical properties of this system [9] . For p ^ 1 
the force is harmonic, F{p) ^ —p while in the oppo- 
site fimit, p ^ 1, F{p) ^ — 1/p- The effective poten- 
tial V{p) = — J^ F{p)dp — (1/2) ln(l -I- p^) is symmetric 
V{p) — V{~p) and V{p) ~ ln(p) when p ^ 1 (we use 
a dimensionless representation for p (3|). The minima of 
the effective potential V{p) is at p = which is of course 
the ideal cooling limit which is not maintained due to the 
fluctuations described here by D. 
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FIG. 1: We show that f{z) = t^'^^°'>W{p,t) {z = p/t^'^) ob- 
tained from numerical solutions of the Fokker-Planck equation 
with a logarithmic potential, converges towards the infinite in- 
variant density [solid blue t = oo curve Eq. (|12|) ]. Thus the 
asymptotic behavior of the system is not globally determined 
by the standard steady state solution. Notice that for t -^ oo, 
f(z) ~ z~^'^ for small z and hence in in the long time limit 
the infinite invariant density f{z) is non-normalizable. Here 
D = 1/2 and initially W{p,t = 0) is a delta function centered 
on the origin. 



Steady state. The steady state solution of W{p,t) is 
found in the usual way: imposing dW steady I dt = from 
Eqs. (|ll2p we have Wsteady ex exp[— y(p)/_D]. This solu- 
tion is normalizable only if £) < 1 and for that case one 
finds the Tsallis distribution [3|, |9[ 



Wsteady {P)=M {I +p') 



2\-l/{2D) 



(3) 



where -^ — ^ (ji)) /[v^r {^^^)] is a normalization con- 
stant. This steady state solution was observed in opti- 
cal lattice experiments [9| where it was shown that this 
behavior is tunable, namely one may control D to ob- 
tain different steady state solutions. Notice that Eq. ([3]) 
exhibits a power law decay for large p which is clearly 
related to the logarithmic potential under investigation. 
From Eq. (jSj we have 



(p'] 



1-3D 



steady 



< £> < 1/3 
l/3< D < 1. 



(4) 



As mentioned in the introduction the behavior 
ip'^) steady = OO is uuphysical since it implies an averaged 
kinetic energy which is infinite jlOJ . 

Bounds on (p^). To start our analysis we consider the 
dynamics of (p^). Multiplying the Fokker-Planck Eq. ([1]) 
with p^ and integrating over p we have, after integrating 
by parts and using the natural boundary condition that 
W{p, t) and its derivative at p — >■ ±oo are zero. 
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where 



/_ ...W{p,t)dp. Obviously we have 



0< 



'l+p2 



<1, 



(6) 



hence 2L» - 2 < 



dip' 

dt 



< 2D, and therefore if we start 



with a finite (p^) we have in the long time limit 
{2D - 2)t < (p2) < 2Dt. 



(7) 



The upper bound clearly implies that (jP) increases at 
most linearly as diffusion persists. The lower bound is 
useful when D > 1 since then it shows that (p^) oc t. We 
call the case D > 1 the diffusive regime. We now turn to 
analyze the cases D > 1 and D < 1 separately since they 
exhibit very different behaviors. 

The case D < 1. We now consider the more interesting 
case D < 1 where a normalizable steady state Eq. (|3]) 
exists. In the long time limit the latter describes well 
the central part of W{p, t) but not its tails which govern 
the growth of (p^) when 1/3 < D < 1 (for Z? < 1/3 
higher order moments diverge and the essential problem 
remains). We employ the scaling ansatz [11| 



wip,t)^ef{p/Vt) 



(8) 



which holds for large p and long t and the exponent a 
will be soon determined. Let us introduce the scaling 
variable z — pjt^l"^ . This is the typical scaling of Brow- 
nian motion, which indicates that for large p diffusion is 
in control; however as we now show, /(z) in far from a 
Gaussian so the process is clearly not simple diffusion. 
Inserting Eq. ([8]) in the Fokker-Planck equation ([1]) and 
using p 3> 1 we find 
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dz^ 



1 



dz 
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/ = 0. 



(9) 



For small z we get / ^ z^^/^ or / ^ z; the latter is 
rejected since f{z) cannot increase with z. To find a 
we require that the small z solution matches the steady 
state, since the latter describes well the density in the 
center. Using Eq. ^ with / ex z~^^^ we have W{p, t) ex 
-l-a+i/{2D)p-i/D ^ffY^Ych is to be matched with the steady 
state solution Eq. ^ Wsteady oc p^^l^ . Hence we have 



1 
2D' 



(10) 



In this case one solution of Eq. (j9]) is immediate: /(z) — 
Az^^'^ . While this solution has the correct small z be- 
havior it does not decay quickly enough at large z |l2|. 
so we need the second solution. The later is found by the 
method of reduction of order: 



/■oo 

f{z) = Az-^"" \ si/^e-'^'/^^ds 

J z 



(11) 
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The variance (p^) . Even though the solution Eq. p2)) 
is non-normahzable, it can be used to find the second 
moment (p^). To see this we introduce a cutoff Pc 
above which our solution Eq. ([T^ is valid. The vari- 
ance is calculated in the usual way exploiting symmetry 



{p')=2 I p^W{p,t)dp . o / .2 



p^W{p,t)dp. (14) 



The first term in Eq. ([H]) is a constant and can be ne- 
glected once the second term is shown to increase with 
time. Inserting the infinite invariant solution Eq. (J12p in 
the second term of Eq. ^1)1 



FIG. 2: In the long time limit the variance (p^) exhibits 
three behaviors: normal diffusion if D > 1, anomalous dif- 
fusion (p2) ~ i3/2-i/(2D) a i/s < D < 1 and steady state 
behavior (p^) — >■ const for D < 1/3. For very short times (not 
shown) (p^) ~ 2Dt since then the influence of the force field 
is negligible (initially particles are on the origin). In the inset 
we show that the ratio of {p^) obtained from simulations to 
(p^) found from the asymptotic formula Eq. (|17p approaches 
unity as time increases. 



The constant A is found by matching the small z solution 
Eq. pT|) to the steady state solution Eq. ([3]). Solving 
the integral in Eq. (|lip we reach our first main result 



f(z) 



Nz 



-l/D /i^jj p 



r(^] 



2D 'AD 



77^ (12) 



where T{a,x) = J^ e^'^s'^^^ds is the incomplete Gamma 
function J13j and T{a) is the Gamma function. For small 
and large z we find 
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Eq. ([T^ is non-normalizable since according to Eq. P^ 
f{z) ~ z~^^^ and hence f^ f{z)'iz = oo. 

Infinite invariant density. We call the non- 

normalizable solution Eq. (112^ an infinite invariant den- 
sity. In Fig. [1] comparison is made between our analytical 
solution Eq. (|12p and numerical solutions of the Fokker- 
Planck equation. As time increases, the solution in the 
scaled coordinate approaches the infinite invariant den- 
sity Eq. ([T^. which describes the asymptotic scaling so- 
lution of the probability density. For any finite long time 
t, expected deviations from the infinite invariant solution 
are found for small values of z (see Fig. [T]). These de- 
viations become negligible at i — >■ oo however they are 
important since they indicate that the pathological di- 
vergence of /(z) on the origin is slowly approached but 
never actually reached, namely the solution is of course 
normalizable for finite measurement times. 
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z^f{z)dz 
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for 1/3 < D < 1. The lower limit in the integral Pc/t^^"^ 
goes to zero when t ^ oo and the diffusion is anomalous 



(p2) ^ 2ti-^ / z'^f{z)dz 



(16) 



Thus the infinite invariant density yields the anomalous 
diffusion in this model. While f{z) ^ z^^/^ for small z 
and is hence non-normalizable, the integral in Eq. (|16p 
is finite: the z^ cures the pathology of the density at the 
origin. Solving the integral in Eq. ([T6l) as well as the 
diffusive regime soon to be discussed we obtain 
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(17) 
For D < 1/3, (p^) is time independent and is determined 
by the steady state solution Eq. 1^. For the intermediate 
regime 1/3 < D < 1 the diffusion is anomalous, while for 
D larger than 1 it is normal in agreement with the simple 
bounds we have found, Eq. ([T]). In Fig. [21 numerical 
solutions for (p^) versus time exhibit convergence towards 
these types of behavior. 

A simple argument for the anomalous diffusion expo- 
nent in Eq. PT)) . for 1/3 < I? < 1, is found by noticing 
that the steady state solution Eq. ([3]) describes the cen- 
ter part of the PDF however with cutoffs determined by 
diffusion (i.e. \p\ < vDt) hence 



ip' 



/Dt 



-VDt 



_P^Wsteadv{p)dp CX 2 / p^-l/D^^ ^ ^ 



Vt 



(18) 

To characterize the distribution of p and to find (p^) ex- 
actly, we need the infinite invariant density which cannot 
be obtained by similar simple scaling arguments. 



The case D > 1. We now investigate the diffusive 
regime D > 1 where the steady state solution is non- 
normahzable. In this case our bound Eq. ([7]) yields dif- 
fusive behavior p (x y/i which suggests a scaling solution 
W{p,t) = t°'f{z) where now a = —1/2 and as before 
z ~ p/^/t 11]. Eq. ^ is still valid and its solution when 
a = -1/2 is 



/W = 



-l/D 



e^nr 



(19) 



(4£))i-'3r(i-/?) 

where /3 = (1 + D)/{2D). No tice th at when D > 1 
we have f{z) ^ exp(— z^/4Z?)/-\/47r_Di which is expected 
from a purely Gaussian diffusive process. Roughly speak- 
ing the potential is responsible for accumulating of par- 
ticles close to the origin which yields the z~'^/^ factor in 
Eq. (|19l) . Similar to the case D < \ the solution Eq. 
P^ exhibits a divergence on z = since f{z) ~ z~^/^. 
However, since now D > 1 the solution is normalizable 
and in this regime we do not find an infinite invariant 
density. Finally with Eq. (IT9l) we obtain {p^) for Z? > 1 
which was given in Eq. (J17p . 

The role of infinite invariant density in physics is now 
briefly discussed. Mathematicians have investigated infi- 
nite invariant measures in the context of ergodic theory 
for many years |14| |. More recently, an infinite invariant 



density 'l5il was shown to describe weak chaos in the well 
known intermittent Pomeau-Manneville map (by weak 
chaos we mean dynamical systems with zero Lyapunov 
exponent which still exhibit stochastic behavior). In the 
context of a model of an electron glass, the distribution of 
eigenvalues of a relaxation matrix were found to be non- 
normalizable which yields slow relaxations and aging [16| . 
Thus three routes to infinite invariant densities are: weak 
chaos, disorder and as we showed here the widely appli- 
cable case of diffusion in a logarithmic potential. These 
vastly different mechanisms all exhibit anomalous diffu- 
sion and the usual ergodic hypothesis breaks down. It 
therefore seems likely and certainly worthy of further in- 
vestigation that a broad ran ge of p hysical systems which 
exhibit anomalous kinetics [17H19| are described by an 
infinite invariant densities. 

Thermal systems. So far we have used the example 
of the motion of two level atoms in an optical lattice, a 
system which is not thermal and its equilibrium distribu- 
tion is strictly speaking not a Boltzmann equilibrium. As 
noted in the introduction we may consider over-damped 
Brownian particles coupled to a thermal heat bath with 
temperature T and get the same results. More precisely 
consider over-damped Brownian motion in the potential 



eo In (a^ 



/2 and diffusion constant D (units m / s). 



With the fluctuation dissipation theorem [20| we have 
dP{x,i) UT d^ d 



dt 



, az.2 + fl™ 1 I z,2^ y^^^l 



eo 



dx"^ dx 1 



(20) 



which after obvious change of notation is the same as Eq. 
dD). In Eq. (pn|) dimensionless time i= e^Dt/a^kbT and 



space X = x/a are used. More importantly our results 
are not limited to one dimension. Indeed an infinite wire 
of radius a, with uniform charge density per unit length A 
yields the logarithmic potential V{r) — Xln{r) for r > a 
and a > must be finite. Such a potential was consid- 
ered by Manning [6| in the context of ion condensation 
on a long polyelectrolyte. Using cylindrical coordinates 
it is not difficult to show that the corresponding three 
dimensional Fokker Planck equation maps onto a one di- 
mension problem (for coordinate r) which is similar to 
ours. We do note that the limit a — > yields new be- 
haviors beyond what we discussed here (this is related 
to the fact that our potential is finite on p = while 
ln{r) — >■ — c» for r — >■ 0, i.e. when a approaches zero). 

Summary. Steady state solutions are commonly as- 
sumed to describe the long time limit of dynamics of 
many thermal and non-thermal systems. This assump- 
tion is one of the pillars on which statistical mechanics 
is built. Therefore it was rewarding to find that for the 
widely applicable process of Brownian motion in a loga- 
rithmic potential, an infinite invariant density describes 
the scaling solution. Both the well known steady state 
solution Eq. ([3]) and the infinite invariant density Eq. 
P^ found here are needed to characterize the long time 
solution. Thus Boltzmann's equilibrium concepts while 
important are clearly not sufficient. 
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